The Glassy Potts Model 
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We introduce a Potts model with quenched, frustrated disorder, that enjoys of a gauge symmetry 
that forbids spontaneous magnetization, and allows the glassy phase to extend from T c down to 
T = 0. We study numerical the 4 dimensional model with q — 4 states. We show the existence of a 
glassy phase, and we characterize it by studying the probability distributions of an order parameter, 
the binder cumulant and the divergence of the overlap susceptibility. We show that the dynamical 
behavior of the system is characterized by aging. 
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The generalization of the Ising model to a frustrated 
model containing quenched disorder, the Ising spin glass, 
has provided us with a large amount of new physics 
0. Replica Symmetry Breaking has been found in the 
Mean Field theory [0 , and mainly numerical simulations 
strongly hint to its validity in finite dimensional disor- 
dered Ising spin glasses 0. 

The need of a generalization of such systems to the 
Potts models has been clear very soon: technical motiva- 
tions are obvious, while physical motivations include the 
need of describing systems where the Zi symmetry of the 
Ising model is not relevant (real glasses being potentially 
among them 0). The most straightforward construction 
of a Potts spin glass, where the spin variables can be in 
q states and are randomly connected by a positive or a 
negative coupling, has been analyzed in detail ||- [fl3|| , 
but it has the unappealing feature (that we will justify in 
the following) of acquiring a spontaneous magnetization 
(i.e. of entering a phase with usual ferromagnetic order- 
ing) at low T values. The glassy regime is only present 
in a small T region, making it unpractical to be studied 
numerically and unplausible for a faithfully description 
of real glasses (that do not order at low temperatures) . 

Here we will define and study what we consider to be 
the "naturally glassy" generalization of the Ising spin 
glasses to the Potts model. We will study the finite di- 
mensional version of the model, and we will show that 
these systems do indeed undergo a phase transition that 
leads them to a glassy phase, different from the usual 
Sherrington Kirkpatrick spin glass phase. 

We regard here as crucial the exact gauge invariance 
that is found in the usual Ising spin glasses (both in 



finite number of dimensions and in the mean field ap- 
proximation). The Hamiltonian of the Ising spin glass 
is — ^ (JiJijUj, where the sum runs over couples of first 
neighboring sites of the (simple hyper-cubic) d dimen- 
sional lattice, the J are ±1 with uniform probability (or 
Gaussian variables), and the spin a take the values ±1. 
Let us consider the site i and transform ct,; — ► — (Tj. If at 
the same time we flip all the 2d couplings Jjj involving 
the site i (and one of its first neighboring sites j) the 
energy of the system does not change. Because of this 
symmetry the expectation value of the magnetization is 
zero, and a ferromagnetic phase is not allowed. 

The generalization to a disordered model of 0- [|l3| is 
not protected by such a local symmetry, and is allowed to 
magnetize (as it indeed does) . We propose instead a gen- 
eralization to a quenched, frustrated spin model, where 
the gauge invariance is preserved. The Hamiltonian of 
our model is 
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where the sum runs over first neighboring sites on a sim- 
ple cubic hyper-lattice (or over all site couples in the 
mean field model), the spin variables a can take q val- 
ues (0, 1, . . ., q — 1), and the Hij are link attached, 
quenched, random permutations of (0, . . . , q — 1) (there 
are ql of them). In the following we will denote this 
model by A4 q \. It is clear that, when written as a sum 
of delta functions, the Ising spin glass has exactly this 
form: in this sense this is a very natural generalization of 
the model, where we just increase the number of allowed 
states. The link (i, j) will give a non-zero contribution to 
the action not, as in the usual Potts model, if a{ = cr,-, but 
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if Gi is equal to IT .j{<Jj). In this model the same gauge 
symmetry we have described before protects us against 
magnetizing: if we transform <Xj from to 1 we will inter- 
change in the quenched random permutations involving 
the site i the state and the state 1. This feature makes 
this model (the glassy Potts model) a good candidate to 
the description of the glassy state. 

The model we have just defined, A4 q \, has a drawback: 
it is very difficult to check if it has reached thermal equi- 
librium. The Zi symmetry of the Ising model is indeed 
precious at this end: in a spin glass checking the symme- 
try of the probability distribution of the overlap is crucial 
for establishing thermalization ||. Our first numerical 
simulations of the model M. q \ have confirmed how dif- 
ficult it is to establish on firm grounds thermalization 
without being able to count on a slow mode that has to 
exhibit a symmetry. It is possible to solve this problem 
at least for even values of the number of allowed states, 
q. One considers a permutation R such that R 2 = 1 
(for example we can change the state 2k with the state 
2k + 1 for k = 0, ... | — 1), and allows in the 6 func- 
tion only permutations that commute with R. We have 
introduced the model A4 co (where co stands for commu- 
tative) where the Hamiltonian of equation (Q) contains 
only permutation II that commute with the permutation 
R = (0, 1, 2, 3) -> (2, 3, 0, 1), i.e. UR = RIL. This model 
is symmetric under R, and invariance under R can be 
tested in order to check if thermal equilibrium has been 
reached. In order to do that we have defined a modified 
overlap w, that is one if two spins are in the same state, 
— 1 for the couples (0,2), (1,3), (2,0) and (3, 1), and zero 
otherwise (q will be the usual overlap, where we sum one 
if two spin are equals and zero otherwise). Because of 
the symmetry we have introduced by selecting only R 
commuting permutations II the probability distribution 
P{oj) is symmetric at equilibrium under lu — > — u). The 
two models are expected and turn out to be equivalent, 
as we will show in the following. When using A4 co it is 
easy to check thermalization, and the coincidence of the 
results with the ones obtained when studying M. q \ shed 
light on their physical meaning. 

We have studied M. q \ and M. C o with q = 4 states in 
4 spatial dimensions d. We have used a normal Monte 
Carlo method. We will present thermalized data in the 
broken phase for lattices of volume L 4 = 4 4 and 5 4 , and 
data in the warm phase for a 8 4 lattice. The data on 
the two smaller lattice volumes have been obtained from 
a slow annealing with ten million full sweeps of the lat- 
tice at each temperature point, the one with L — 8 have 
used one million sweeps per T point. We have averaged 
over ten realizations of the disorder. We have preferred 
to have long thermal runs, since thermalization of the 
samples needs to be completely sure to make the results 
reliable, and to keep the number of samples quite small. 
The numerical simulations have taken of the order of two 
months of medium size workstations. We only report re- 



sults for which we are sure of having reached full ther- 
malization (the main criterion used being the symmetry 
of the probability distribution of M co , and the request 
of a good stability in time of the observables). 

Working on lattices of linear size 4 and 5 we have suc- 
ceeded to get some control over the finite size behavior. 
A large amount of evidence, that we will describe in the 
following, makes clear the existence of a phase transition 
to a low T glassy phase. 

The Binder parameters of the modified overlap 
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shows the clear signature of a phase transition. For M. co 
~ with good accuracy in the warm phase (because 
of the symmetry we have implemented) . It becomes dif- 
ferent from zero at T ~ 1.4, and grows basically linearly 
(in our statistical accuracy) at low T. At T — 1.2 (the 
lowest T value we are sure we have thermalized both at 
L = 4 and L = 5) g u ~ \. In the precision given by our 
statistical errors (of the average over disorder), evaluated 
directly for the Binder parameter by a jack-knife method, 
the L — 4 and L — 5 results coincide (at T ~ 1.2 the rel- 
ative error on the Binder parameter is of he order of ten 
percent: for example for L — 5 we have g^ = 0.45±0.10). 
For the model M. q \ g u has the same pattern but it be- 
comes slightly negative at T ~ 1.7 (but close to zero): at 
low T it increases like for the other model. The Binder 
parameter of the straightforward overlap does not behave 
in a interesting way for both models. From the analysis 
of the Binder parameters we deduce as a first guess that 
T c ~ 1.5. 

In figure (|l|) we plot the Replica Symmetry Breaking 
parameter introduced in MM, 
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for M.co (where we have assumed that (u>) — 0. p seems 
to give an even clearer signature of the phase transition 
than the Binder parameter (qualifying it without ambi- 
guities as a Replica Symmetry Breaking transition). In 
the infinite volume limit p is zero if the overlap distri- 
bution is self-averaging, and becomes non-zero if Broken 
Replica Symmetry makes it non-self-averaging. Refer- 
ence jl4| shows that while the Binder parameter is very 
effective in detecting symmetry breaking accompanied by 
breaking of spin reversal symmetry, but when spin rever- 
sal symmetry is absent p tends to be a better estimator. 
Figure (]!]) shows a sharp transition (notice that the ver- 
tical scale is logarithmic). Again a good estimate for T c 
is 1.5. The evidence for a phase transition is clear. The 
fact we can detect it by using p makes clear it is a Replica 
Symmetry Breaking transition. Again, the situation in 
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M q \ is very similar: there is a change of regime close to 
T = 1.5, where p does not go to zero with increasing L. 

Maybe the most interesting evidence about the behav- 
ior of the system comes from the analysis of the proba- 
bility distributions of the overlap averaged over the dif- 
ferent disorder samples. At high T P(lu) in M. q \ tends to 
a Gaussian when increasing the lattice size. We show in 
figure @ P{u) at T = 1.25 for L = 4 and 5. There is a 
clear non-trivial behavior (notice that the support in lu is 
very extended, i.e. P(ui) is very flat). On larger lattices a 
peak in lu ~ is emerging, separated by a flat minimum 
from a second peak (the minimum becomes sharper at 
lower T values, but there we know we did not reach full 
thermal equilibrium and we cannot safely attach a firm 
significance to the data): this behavior is reminiscent of 
the one of REM models, and fits with what we expect for 
a glassy state The important point here is that we 
have been able to thermalize the L — 5 system even in 
the cold phase (T = 1.25 is the lowest temperature value 
where we are sure about an adequate thermalization) . 

It is useful (and necessary) to look at the individual, 
single sample probability distributions P j{oS) in order to 
qualify the behavior of the individual systems. In fig- 
ure (||) we plot Nj(u>) (the non-normalized Pj{lu)) for 
four typical samples, versus u> at T = 1.25. The level 
of the asymmetry of the histograms is a measure of our 
statistical error (and the fact the functions look sym- 
metric a sign of a good thermalization). The Nj(lu) are 
non-trivial: some samples have double peaks, some have 
their support close to lo = 0, some have support at zero 
overlap and peaks at finite overlap. On a qualitative level 
we remark that system looks harder than the usual spin 
glass: the dynamics is more jumpy, visiting in a quite 
discontinuous manner different parts of the phase space 
(that is why checking thermalization has been difficult 
and crucial). Our evidence seems to suggest that the 
free energy phase space is golf course like: deep minima 
do not have large basins of attraction. 

The probability distribution for the modified overlap 
lu in the model A4 q \ is similar to the one we have shown. 
The detailed shape is not exactly the same, but it also 
becomes non-trivial at T — 1.25. On the L = 5 lattice 
a two peak structure starts to emerge from T = 1.25 
down (one in lu ~ and one at a finite lu value). The 
non-modified overlap distributions for both models enjoy 
the same main features (even if they are non-symmetric, 
and, as we have already discussed, the thermalization is 
better checked by using the symmetry of P{lu) in M. co )'- 
it tends to a Gaussian when increasing L at T = 1/0.6, 
it is non-trivial at T — 1/0.7, and it develops a two peak 
structure at lower T values. 

We have used the T data in the warm phase from the 
large, L = 8 lattice to fit the divergence of the overlap 
susceptibility x w = V(tu 2 ), that T — > T+ behaves as 
(T — T c ) 1 . From our data we can only give a preliminary 
estimate, that puts T c among 1.4 and 1.5 (compatible 



with the value we have deduced from the direct analysis 
of P(lu)) and the exponent 7 among 1.3 and 1.5. This 
value is different from the one quoted for the Ising 4D 
spin glass, 7 = 2.10 jl5). 

The dynamical behavior of the system shows all the 
typical features of the complex dynamics. In figure (^) 
we show the aging behavior: the spin-spin correlation 
function C{t,t w ), depending on the waiting time t w and 
on the time t, versus t for different waiting times t w . The 
rate of the time decay depends on t w . 

Let us notice at last that the energies of the two models 
are very similar: on the L = 5 lattice they are equal in 
our statistical precision, while on the L = 8 lattice they 
are of the order of one per one thousand. The two models 
seem to present the same kind of critical (and off-critical) 
behavior, and our best guess is that they do belong to 
the same universality class. 

We have introduced a disordered generalization of the 
Potts model that we regard as a very hopefully candidate 
to describe the glassy state. Our numerical simulations 
of the 4 state, 4 dimensional model show clearly the exis- 
tence of a glassy phase, and they stress large differences 
with the usual Sherrington Kirkpatrick spin glass phase, 
exhibiting a more discontinuous behavior, reminiscent of 
the Random Energy Model one step Replica Symmetry 
Breaking Q. 

We thank Renate Loll for the gift of a nifty script com- 
puting commuting permutations. 
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FIG. 2. P(u>) (symmetrized) versus uj at T = 1.25. An- 
nealed runs with ten million sweeps per T point. Solid line 
for L = 4 and dashed line for L — 5. 
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FIG. 3. Nj(lu) (the non-normalized Pj(ui)) for four typ- 
ical samples, versus u) at T = 1.25. Annealed runs with ten 
million sweeps per T point. 



FIG. 1. p(T) versus T. Solid line for L — 4 and dashed line 
for L — 5. Notice the vertical logarithmic scale. 
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FIG. 4. Dynamical correlation function C(t,t m ) versus t 
for different values of t w . 
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